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Abstract. We prove universality at the edge of the spectrum for unitary (/9 = 
2), orthogonal (/3 = 1) and symplectic (/3 = 4) ensembles of random matrices in 
the scaling limit for a class of weights w{x) = e~^(^' where V is a polynomial, 
V{x) = K2m2:^™ + ■ ■ • , K2m > 0. The precise statement of our results is given 
in Theorem 11.11 and Corollaries 11.21 ll.iji below. For a proof of universality in 
the bulk of the spectrum, for the same class of weights, for unitary ensembles 
see DKMVZ2 , and for orthogonal and symplectic ensembles see DG . 

Our starting point in the unitary case is [pKMVZ2 , and for the orthogonal 
and symplectic cases we rely on our recent work .DG. . which in turn depends on 
the earlier work of Widom ^ and Tracy and Widom ITW2I . As in IDG) , the 
uniform Plancherel— Rotach type asymptotics for the orthogonal polynomials 
found in )DKMVZ2I plays a central role. 

The formulae in IWI express the correlation kernels for /3 = 1 and 4 as a 
sum of a Christoffel-Darboux (CD) term, as in the case (3 = 2, together with a 
correction term. In the bulk scaling limit IDGI . the correction term is of lower 
order and does not contribute to the limiting form of the correlation kernel. 
By contrast, in the edge scaling limit considered here, the CD term and the 
correction term contribute to the same order: this leads to additional technical 
difficulties over and above IDGI . 



1. Introduction 

This paper is a continuation of |DGj . In jPGj . the authors proved universahty 
in the bulk for orthogonal and symplectic ensembles: here we prove universality at 
the edge for orthogonal and symplectic ensembles, and also for unitary ensembles. 
For the convenience of the reader, and to fix notation, we now summarize some 
of the basic theory of invariant ensembles (/3 = 1, 2 or 4), borrowing freely and 
extensively from the introduction in DG]. We are concerned with ensembles of 
matrices {A/} with probability distributions 

(1.1) Vn p{M)dM = -^e-''^^^^UM, 

Zn,I3 

for f] = 1, 2 and 4, the so-called Orthogonal, Unitary and Symplectic ensembles, 
respectively (see |M1| V For /? = 1, 2, 4, the ensemble consists of N x N real 
symmetric matrices, N x N Hermitian matrices, and 2N x 2N Hermitian self- 
dual matrices, respectively. In general the potential Vf}{x) is a real- valued function 
growing sufficiently rapidly as \x\ oo, but we will restrict our attention henceforth 
to Vg's which are polynomials. 



(1.2) 



Vl3{x) = K2m,pX^"' H , K2m,/3 > 0. 
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In , dM denotes Lebesgue measure on the algebraically independent entries of 
M, and 2m^i3 is a normalization constant. The above terminology for /3 = 1, 2 and 
4 reflects the fact that Hl.l|l is invariant under conjugation of M, M i— > UMU^^, 
by orthogonal, unitary and unitary-symplectic matrices U. It follows from (|l.lfl 
that the distribution of the eigenvalues xi, • • • , x^r of Af is given (see IMll ) by 

1 ^ 

(1.3) PN.fsixi,- ■ ■ ,Xn) ^ ^ — TT \Xj - XkfYl'^Pi^j) 

'P l<j<k<N j=l 

where again ZN,f3 is a normalization constant (partition function). Here 

''e-^'3(^), (3=1,2 



(1.4) W[3{x) 



-2V,(x)^ /J = 4. 



(The factor 2 in z«^=4 reflects the fact that the eigenvalues of self-dual Hermitian 
matrices come in pairs.) Let {pj}j>o be the normalized orthogonal polynomials 
(OP's) on R with respect to the weight w = w^=2, and define = pjW^^"^. Note 
that {(l)j,cj)k) = Sjk where (•, •) denotes the standard inner product in L^(IR). 

For the unitary matrix ensembles an important role is played by the Christoffel- 
Darboux (CD) kernel 

(1-5) KN{x,y) = KN,2{x,y) = ^ (l)k{x)(l)k{y)- 

k=0 

In particular the probability density (|1.3|l . the ^-point correlation function Rn,i,2 
and also the gap probability i?2(0; J) that a set J contains no eigenvalues, can all 
be expressed in terms of Kn, see e.g. |M1| . For example 



(1.6) Rnj,2{xi,--- ,xi) = det{KN{xj,Xk))i<j,k<i- 

The Universality Conjecture, in our situation, states that the limiting statistical 
behavior of the eigenvalues xi, • • • ,xjv distributed according to the law H1.3(l . in 
the appropriate scale as ^ oo, should be independent of the weight w^, and 
should depend only on the invariance properties of 'PN,f3, /3 = 1, 2 or 4, mentioned 
above. Universality has been considered extensively in the physics literature, see 
e.g. lEFTllHel lHWe SeVe . 

The kernel Kn{x, y) can also be expressed via the Christoffel-Darboux formula 

(1.7) KN{x,y)=bN-i , 

x-y 

where bN-i is a coefhcient in the three-term recurrence relation for OP's, see |Sz|. 
In view of the preceding remarks it follows that in the case /3 = 2, the study of the 
large N behavior of Pn,2, and in particular the proof of universality, reduces to the 
asymptotic analysis of b^-i and the OP's pat+j- with j = or —1. By a fundamental 
observation of Fokas, Its and Kitaev [FoIKi] the OP's solve a Riemann-Hilbert prob- 
lem (RHP) of a type that is amenable to the steepest descent method introduced by 
Deift and Zhou in |DZj and further developed in [PVZ] . In |l)KMVZ1llT)KMV/2j 
the authors analyzed the asymptotics of OP's for very general classes of weights. 
In particular they proved the Universality Conjecture in the bulk in the case (3 = 2 
for weights 'w{x) = e~^^^^ where V{x) is a polynomial as above, and also for 
■w{x) = e^^^(^) where V{x) is real analytic and V{x)l\og \x\ -foo, as oo. 
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The bulk scaling limit as — > oo is described in terms of the so-called sine kernel 
Koo{x - y) where 

„N / N sinTri 

(1.8) i^oo(t) = -—■ 

nt 

For example |DKMV!Z2l Theorem 1.4], for w{x) = e~^'^^\ V{x) polynomial, and 

for any ^ = 2, 3, • • • and r, yi, • • • , in a compact set, one has as iV — » oo 

(1.9) 

W^^kw ^ ^^^y ^^^^ ^ dct(if^(y, - y.))l<.,.•<^ 

The scale a; = y / K]^{Q,0) is chosen so that the expected number of eigenvalues 
per unit y-interval is one. This scaling in the bulk is standard in Random Matrix 
Theory. Indeed for any Borel set B cM., 

(1.10) / Rn,i=i,2{x) dx = E{ number of eigenvalues in B }. 

Jb 

Thus by (|l.t)|) Kjq[Q,Q) — i?jv, 1,2(0) gives the density of the expected number 
of eigenvalues near zero. From (|1.9|l . we see that, in the appropriate scale, the 
large N behavior of the eigenvalues is universal (i.e. independent of V). Pioneering 
mathematical work on the Universality Conjecture in the bulk was done in |PS| and 
for the case of quartic two-interval potential V{x) = N{x'^ — tx'^), t > (sufficiently) 
large, in |^. We note again that all these results apply only in the case /? = 2. 

In the case (3=1 and 4 the situation is more complicated. In place of (|1.5|l one 
must use 2x2 matrix kernels (see e.g. IMIIITW^ I 
(1.11) 

(SN,i{x,y) {SN.iD){x,y) 
I , N even, 
{eSN,i){x,y) - ^sgn{x - y) SN,iiy,x) 

and 

^ I SnA^^v) {SN,iD){x,ij) 



(1.12) KnA^^v) 2 



Here SN.f}{x,y), (3 = 1,4, are certain scalar kernels (see (|1.17ll . (|1.18|) below), 
D denotes the differentiation operator, and e is the operator with kernel e{x,y) = 
^ sgn(a;— y)^. Such matrix kernels were first introduced by Dyson |Dy| in the context 
of circular ensembles with a view to computing correlation functions. Dyson's 
approach was extended to Hermitian ensembles, first by Mehta | M2 | for V{x) = x^, 
and then for more general weights by Mahoux and Mehta in |MaM| . A more direct 
and unifying approach to the results of Dyson-Mahoux-Mehta was given by Tracy 
and Widom in ^TW2j . where formulae p. 17(1 . H1.18|l below were derived. We see 
that once the kernels SN,ii{x^y) are known, then so are the other kernels in K^jj- 
As in the case (3 = 2, the kernels K^^p give rise to explicit formulae for Rn,i.i3 and 
i?/3(0; J). For example for /3 = 1,4 

(1.13) Rn,i,i3{x) =Ri^i3{x) = ]^ixKN,i3{x,x) 



"'^We use the standard notation sgnx = 1, 0, —1 for x > 0, x = 0, x < 0, respectively. 
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and 

RN,2.^p{x,y) = ^ {tiKN,i3{x,x)){tiKN,i3{y,y)) - ^tr {KN,0{x,y)KN,p{y,x)) , 

and so on, see |TW2j . We will discuss some of the hterature on edge scahng after the 
statement of our results, Theorem ll.il Corollary II .21 and II .31 below. As indicated 
above, formula Hl.ll|l only applies to the case when N is even. When N is odd, 
there is a similar, but slightly more complicated, formula (see |AFNvM] ) . As in 
|DG| . throughout this paper, for P = 1, we will restrict our attention to the case 
when N is even. We expect that the methods in this paper also extend to the 
case (3 — I, N odd, and we plan to consider this situation in a later publication. 
Of course, in situations where the asymptotics of Hl.ll|l has been analyzed (e.g. 
V{x) — x'^) for all iV as ^ oo, the limiting behavior of Rn,i,i3=i is indeed seen 
to be independent of the parity of N (see e.g. jMll INW) ). 

Let {qj (a;)}j>o be any sequence of polynomials of exact degree j, qj {x) — qj.jX^ + 
• • • , gfjj ^ 0. For j = 0, 1,2, - • •, set 



(1.14) ^,,0(x) = 



qj{x)wi{x), /3 = 1 

q,{x){w4{x)y/^ (3^4. 



Let Mat^i denote the N x N matrix with entries 

(1.15) {MN,i)jk^ {^j,i,e^ks), 0<j,k<N-l, 
and let Mjv,4 denote the 2A x 2A^ matrix with entries 

(1.16) {MNA]jk = {'^j,i,D^uA), Q<j.k<2N -I, 

where again (•, •) denotes the standard real inner product on R. The matrices M^^i 
and Mn^4 are invertible (see e.g. jAvMl (4.17), (4.20)]). Let hn,i, f^N,4 denote 
the inverses of Mm,i, Mm, a respectively. With these notations we have |TW2| the 
following formulae for Sn,p in p. 11(1 . H1.12|l 

N-l 

(1-17) SNs{x,y) = - '^j,i{x)if^N,i)jk{£ipk,i){y) 

j,k=a 

2N-1 

(1-18) SnA^'V)"^ X! '^jAi^) if^NA)3k^k,4iy)- 

j,k=0 

An essential feature of the above formulae is that the polynomials {qj} are arbitrary 
and we are free to choose them conveniently to facilitate the asymptotic analysis 
of l|l.ll() . (|1.12|l as A^ ^ cx) (see discussion in |DGj and H1.21|l below). 

In order to state our main result we need more notation. For any m € N let 
V{x) be a polynomial of degree 2m 

(1.19) V{x)=K2mX^"' + -- - , K2m>0 

and let w{x) = wp^2{x) — e^^*^^-* as before. Let Pj{x), j > 0, denote the OP's with 
respect to w, and set 4>j{x) = pj{x){w{x)y^'^ , j > 0, as above. For /3 = 1,4 set 

(1.20) Vfsix) EE ^V{x) 
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and let N be even. Then by H1.4|l . Wi = e ^Vi — g v ^^^^ y;^ = g = e 
This ensures that for the choice qj = pj in (|1.14|l 

(1.21) 1pj.p=i{x) = 1pj.f3=4{x) = (l)j{x), 

which enables us in turn to handle S^^i and iS'jv/2,4 in p.l7|l . (|1.18|) simultaneously 
(see |DGI Remark 1.3]). Henceforth and throughout the paper, Kp^ denotes the 
Christoffel-Darboux (CD) kernel (II. 5|) . (|1.7|) constructed out of these functions (pj. 

For the bulk scaling hmit in DKMVZl (/? = 2) and [021 W = 1,4), the authors 
used the standard scale of one (expected) eigenvalue per unit interval. At the edge 
it is standard (see e.g. jXfflSi) to use a slightly different scaling which ensures that 
the kernel i^Aiiy(C, v) (see l)1.25|l below) appears in the limiting forms H1.26|l . H1.27|l . 
(|1.28|l below, without any additional factors. Note that formula H1.10|l also holds 
for /3 = 1, 4 and so RN,i=i.pix) gives the density of the expected number of (simple) 
eigenvalues near x for /3 = 1, 2, 4. In view of (|1.10|) . and also in view of H1.13|l and 

frm . fm 

(1.22) ^ 

Rn,i,2{x) ^Kn{x,x), Rn,iAx) ^SnA^.x), Rn/2.iAx) ^ -Sn/2,4{x,x). 

To leading order, the right edge of the spectrum is located at the point cn + rfjv 
where cn, rfjv are the Mhaskar-Rakhmanov-Saff numbers in H3.1() . (|3.2f) below. For 
all three cases, in the neighborhood of cjv + d,N, we use the scale 

(1.23) ^^^W^g^(i + ^^)+rf^ 



where a at is given in (|3. 10(1 (2') below. As we will see (cf. Remark ll.3l below) this 
scaling differs slightly from a scale of one (expected) eigenvalue per unit interval. 

It turns out that the off-diagonal elements in K^^p scale differently as iV ^ cx). 
On the other hand, the statistics of the ensembles are clearly invariant (cf. discussion 
following H2.8|) below) under the conjugation 

for any scalar A. For example, this is obviously true for the cluster functions T^^i^p, 
/3 = 1 or 4, which have the form 
(1.24) ^ 

a 

where the sum is taken over all permutations of {1, • • • , 1} (see |TW2I p. 816]), etc. 
Denote 



(1.25) 



Set 



^ . Ai(OAi^(r?)-Ai'(OAi(yy) 

Ai(z + ^) Ai(z + 77) dz. 



Theorem ll.il and CoroUarv 11.21 and 11.31 below are the main results in this paper. 
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Theorem 1.1. Let (5^2, 1 or 4. For any V{x) of degree 2m as in H1.19|l define 
Vf3{x) and wp^x) as in (|1.2Q(I . (|1.4|l . Fix a number Lq- Then there exists c — 
c(Lo) > such that as N ^ oo^ the following holds uniformly for ^,77 G [Lq, +00). 
In the case (3 — 2: 

(1.26) En,2 = ^K^i^^^'Kv^''^) - KA^ry{i,r^) ^ 0. 
In the case (3 — \ : 

(1.27) ^ T^if^Y'^(^(^\ryW) - K^\i,v) - 
wh- 



ere 



1 



(i^(l))n(e,^) = {K^'^)22{V,0 ^ KA^ry{^,v) + ^ ^KO ' J Ai(t) 

(i^(i))i2(e,r?) ^ -d.KA^ryitri) - ^ ^{0 Ai{rj) 

KA^ry{t,v)dt 

--j^M{t)dt+-j^ ki{t)dt- j Aiit)dt--sgn{C~v)- 
In the case (3 — 4: 

(1.28) £na ^ (C^'^^^^"^^) - K^'Hi,v) -> 
where 

2(ifW)^^(^^^) ^ 2(ifW)22(,y,0 = ^A..,(e,^) - 2 • m)dt 

2{K^^^),2{tv) = -dr^KArryii,!!) - ^ Ai(e) Ai(,7) 

/OO -1 /'OO nOO 

KA^ry{t,v)dt+- J Ai{t) dt ■ J Ai{t) dt. 

For the error term we have as N —^ oo 

(1.29) £n,i = o(l) 

1 

£■^^,4 = o(l)e'"«e-"' 
uniformly for ^,r] € [Lq, +oo). 

Remark 1.1. For /? = 4, but not for /3 = 1, our methods actually prove that 
Sna = 0{N~'^/'''^™-^)e~'^^e~^^ . In order to obtain power law decay for £n,i, it would 
be sufficient to obtain power law decay in the error term in |D(tI Theorem 2.2]: such 
power law decay can be obtained using more sophisticated estimates as in [DGKVj . 



^For /3 = 1, 4, Af is even. 
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We immediately have the following result. Recall formula H1.24|l for the cluster 
functions for /3 = 1,4; for (3 — 2, the cluster functions have the form |TW2I p. 815] 

TN,lAyi'--- ,yi) = J^KN{ya,,ya2)KN{ya2,ya3) - ■ ■KN{yai,y<7i)- 
a 

Corollary 1.2. Let P = 2, 1 or 4. Let V be a polynomial of degree 2m and let 
K''^\ 13 ^ 1,4: be as in Theorem \L1[ Fix a number Lq . Then for /3 = 1 and 
I — 2,3, ■ ■ ■ we have uniformly for ^i, • • • , i^; > Lq 

hm ^—T^,M^,r\---,i^,rA 

(1.30) ^ ' 



For (3 — A, the same result is true provided we replace TV ; i 1^/2,;, 4 o,nd K^^^ — > 
ifW , For P ^ 2, the same result is true provided we replace Tjv,i,i Tm,i,2, 
K^^^ KAiry, 51 ~* jj '^^'^ remove the trace. 

Together with some additional estimates (see Section|2Il, Theorem ll . II also yields 
the following universality result for the gap probabilities. Recall that for a 2 x 2 
block operator A = (^y)ij=i,2 with ^11,^22 in trace class and Ayi,Ai\ Hilbert- 
Schmidt, the regularized 2-determinant (see e.g. is defined by det2(/ + A) = 
det((/ + A)e-^)e*'-(^ii+^22). 

Let Ai denote the largest eigenvalue of a random matrix M . 

Corollary 1.3. Let P = 2, 1 or 4. Let V be a polynomial of degree 2m and let 
K'-I^\ P^ 1,4 be as in Theorem \L1\ Fix a number Lq . Then the following holds. 
In the case P — 2: 

(1.31) Jim^Prob{Ai < (Lo)^} = det(/ - ^ F'^^\L,). 
In the case P ~ 4: 

(1.32) y^J^oh{\, < (L„)W} = ^det(/-ifW)|^,(j^^^^^^3)) ^ F('\Lo). 



In the case P — 1, let g{£_) — \/4-\-£^ , G ~ diag(g,g ^). Then 
(1.33) 



hm^ProbjAi < (Lo)^^)} = ./det2(/ - G/^(i)G~i|,,,.. ^ F«(io). 



Remark 1.2. The regularized 2-determinant is needed for P = 1 because the opera- 
tor with kernel ^ sgn(^ — 77) is Hilbert-Schmidt but not trace class in L^([Lo, -l-oo)). 
The auxiliary function g is needed to ensure that GK^^^G~^ indeed has a 2- 
determinant: there is considerable freedom in the choice of the function g, see 
Remark 12.21 below. 

Remark 1.3. From Theorem 11.11 and 11.22|1 we have as — *■ 00, 



Cat 



(1.34) ^^-^RN,iA{t^^>) =KMry{t,t) + -Ai{t) J Aiiu)du + o{l) 



— 00 
00 



'''' TRN/2,lA{t^''^) = l^Airy(i, <) - ^ Ai{t) I Ai{u) du + o(l) 
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uniformly for t in any fixed half- line [Lo,+oo). In particular the density of the 
expected number of eigenvalues at the edge of the spectrum cn + (In is given by 

72 = (Ai'(0))2 = 0.066987484 



for the indicated values of /3 = 2, 1, 4, where we have used the formula KAiry{t, t) = 



{Ai'{t) f - t{Ai{t))^ and f_^Ai{u)du = |, f^Ai{u)du = i (see IXbSt] ). Thus 



setting t t/jp, P — 2, 1,4, rescales the axis so that the density of the expected 
number of eigenvalues per unit i-interval is one. 

The distributions F^^^Lq), (3 — 1,2,4, are the celebrated Tracy-Widom distri- 
butions which turn out to have applications in an extraordinary variety of different 
areas of pure and applied mathematics (see for example the recent review |TW6| V 
The distributions F(^)(Lo) can all be expressed in terms of a certain solution of 
the Painleve II equation ( [TWl ITWE] 1 . 

The literature on edge scaling, in particular in the physics community, is vast, 
and we make no attempt to present an exhaustive survey. Rather we will focus on 
aspects of the literature which are particularly relevant to this paper. In the physics 
literature, early work on edge scaling for /? = 2 is due to Moore |Mo| and Bowick 
and Brezin |BoBr| . In the mathematical literature for (3 — 2 with Gaussian weight 
V{x) = x^, early work can be found in Forrester [F| and in the seminal work of 
Tracy and Widom |TW4j , where the authors derived the Painleve II representation 
mentioned above for F^^\ For (3=1 and 4 in the Gaussian case V{x) = x^, the 
Painleve expressions for F^-^'> were obtained by Tracy and Widom in |TW5j . but 
without computing directly the edge scaling limit of the Fredholm determinants. 
The edge scaling limits of matrix kernels K^^p, (3 = 1,4, in the Gaussian case 
were obtained by Forrester, Nagao and Honner in |FNH| . The convergence of the 
Fredholm determinants in the Gaussian case for /3 = 1,4 (and also for (3 = 2) was 
first proved only recently by Tracy and Widom in |T W3j . 

Universality at the edge for (3 = 2 was considered by many authors in the physics 
literature (see e.g. IKaFrj l. and for the cases /? = 1,4 see e.g. |Se Ve| . The proof of 
universality at the edge for = 2 in Theorem 11.11 above is based on the estimates 
in |DKMV!Z2] and does not use any results from jWI ITW2I IDCtI . Many researchers 
have noted that universahty at the edge for /3 = 2 is true (see e.g. |CKu| ). but we 
believe that the details of the proof (Theorem ll.il (3 = 2) have not been written 
down previously. In |StlllSt2IISt3j . for (3 = 2,1, 4, Stojanovic proves universality at 
the edge (and also in the bulk) in the special case of an even quartic (two-interval) 
potential considered previously by Bleher and Its |BI| lor (3 = 2. Stojanovic uses 
a variant of the formulae in |^ together with the asymptotics for CP's obtained 
in BI,. Universality for the distribution of the largest eigenvalue for a wide class 
of real and complex Wigner ensembles (see jMlp was proven by Soshnikov in |So| : 
the methods in |So| are completely different from those in the present paper and 
are based on the method of moments. Laguerre ensembles have been considered by 
many authors, see e.g. |f1 IFNH) . Various universality issues at the soft edge, and 
also at the hard edge and in the bulk, for generalized Laguerre ensembles for /3 = 2 
were analyzed recently in The authors are currently completing an analysis of 



(1.35) 



74 



7i 



(Ai'(0))2 + - Ai(0) ^ 0.185330168 
i(Ai'(0))2 - -L Ai(0) = 0.001954035 
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universality questions for such ensembles in the cases (3 — 1 and 4, together with 
Kriecherbauer and Vanlessen, see jPGKVj . 

We complete this introduction with a description of Widom's result fWl which 
is basic for our approach in this paper. Widom's method applies to general weights 
Wfj with the property that w'^/wp is a rational function. This property certainly 
holds for our weights as in H1.4|l . (|1.2(l . and also for general Laguerre type weights 
which we consider in the forthcoming paper |DGKVj . Introduce the matrices 

(1.36) Dn = {{D(l)j,(j)k))o<],k<N-i, eN = {{e(t)j,4>k))o<],k<N-i- 

It follows from jTWll Section 6] that the matrix Djv is banded with bandwidth 
2n + 1 where 

(1.37) ?i=2to-1. 

Thus (Z?Ar)jfc = if \j ~ k\ > n. Next, let N be greater than n, and introduce the 
following A^-dependent n-column vectors 

$i(x) = (0Ar_„(a;),- 

$2(2;) = (0Ar(x),- • • 

e^i{x) = (e(j)N-n{x), 
e^2{x) = (e07v(a;), • • 
and the following 2n x 2n matrices consisting of four n x n blocks 

{{e4>j , 4>k)) N -n<j.k<N +n-l- 



(1.38) 



,e(f>N+n-lix))'^ 



(1.39) 
and 
(1.40) 

where 



B 




Bo 



A = 






{{D(j)j , (j)k))N-n<j,k<N+n-i- Finally, set 

/„ + {BA)n 
{BA)2i 





C12 ^ 




V C21 


C22 1 


H 




Note that 
(1.41) 



Cii 



„ + B12A21 — In — Bi2D21- 

The main result in W is the following pair of formulae for Sn,i and S^/2A 
SnAx, y) = Kn(x, y) - ($i(x)^, 0^) • {ACihn - BACy^f 

■{e^,iyf,e^2iyff 



(1.42) 
and 
(1.43) 



5'Ar/2,4(a;,y) = K^ix^y) + <i>2ixf ■ D21 ■ e^i{y) 

+ <^>2ixf ■ D2iC^^BiiDi2 ■ £$2(2/). 
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Observe that Sn,i and S'jv/2,4 are sums of the [3 — 2 kernel KM{x^y) together with 
correction terms that depend only on (j^N+j for j G {— • • • — 1}. The /3 = 4 
case is different from the case (3—1 since, by (|1.18|l . for any a; G M, 

(1.44) S'Ar/2,4(2^,+oo) 0, KN{x,+(yo) = Q. 
Therefore in H1.43|l for any (even) N and for all a; G M 

(1.45) $2(2:)^ • £"21 • e$i(+oo) + <^2{xf ■ D2iC{iBnDi2 ■ e$2(+oo) = 0. 

As the entries of $2(2^) are functionally independent, and as D12 is invertiblc for 
large N (see 1021 (2.13)]), it follows that 

(1.46) e$i(+oo) + C^i^BuDu ■ e«'2(+oo) = 
for large N . From the definition of e for any integrable ip 

(1.47) eij{y)^- tp{t) dt - l ^/j{t) dt ^ €i/j{+oo) - tp{t)dt. 



2 . 

Hence (|1.43|) . (|1.45|l imply 

<^i{t)dt 



(1.48) 



5*^/2,4(2;, y) = Kn{x, y) + $2(2;)^ • L'21 • ^ 



y 



+ ^2{xf ■ D2iC^,^BnDi2 • ( - / *2(i) dt^ . 



Formula p.48|l makes clear the decay properties of Sn/2a{x, y) a,s x,y ^ +00. Note 
that Sn,i does not satisfy (|1.44|) : this is the reason why we introduce auxiliary 
functions (cf. G — dia,g(g, g^^)) when proving convergence of the determinant in 
Corollarv ll.3l As noted earlier, the question of convergence of the determinants for 
/? = 1, 4 in the Gaussian case was first treated in |TW3| . 

The following observations apply to the 21 entries in the matrix kernels in the 
P —I and 4 cases. Note that by H1.17|l . {eSN^i){x,y) is skew symmetric. Thus 

(1.49) {£SN,i){x,y) ^ {£SN,i){x,y) ~ {eSN,i){y,y) ^ ~ [ SN,i{t,y)dt. 



Also, from (|1.18|l . we see that {eSN/2A){+oo,y) = for all y G M. Together with 
(|1.47|l . this implies that 



(1.50) ieSN/2A){x,y) = - SM/2A{t,y) dt. 

J X 

These observations simplify evaluation of integrals of the CD kernel, and also inte- 
grals of the functions (f)N+j in Sections Eland 01 below. 

Remark 1.4. We note that (|1.49|) is also true for S'Ar/2,4, but p.50|l is more relevant 
for the calculations that follow. 

In Sectional we prove Theorem 11.11 and CoroUarv II . 31 using results on the edge 
scaling limits of the CD terms and the correction terms in Kjq^i and Kiq^^. These 
scaling limits are proved in turn in Section |3| for the CD terms, and in Section^ 
for the correction terms. Note that CoroUarv 1 1.21 is an immediate consequence of 
Theorem ll.il 

Notational remark: Throughout this paper c, c', C, C(m), ci, C2, • • • refer to 
constants independent of N, ^, 77. The symbols c, c', C, • • • refer to generic constants, 
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whose precise value may change from one inequahty to another. The symbol cn 
however always refers to the iV-dependent constant below. 
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2. Proofs of Theorem II. II and Corollary ^31 

The key estimates for the proofs of Theorem 11.11 and CoroUarv 1 1 . 31 are obtained 
below in Section for the CD terms and in Section 0] for the correction terms. 



2.1. Proof of Theorem ll.il Inequality (|3.8fl proves the result for the (3 = 2 case. 

In the case /3 = 4, we use (|1.48|) and consider the CD part and the correction term 
separately. The properly scaled 11, 22 and 12 entries of Kj^^^]^ V^^'') converge 

to the corresponding entries in H1.28|l et seq. with the error estimate o{l)e~'^^e~'''^ , 
uniformly for ^,7/ G [Lq, +oo): this follows from (|3.8() for the CD kernel part, and 
from (|4.22|l and H4.17(l , respectively, for the correction term. By H1.50|l , 1)3.56(1 and 
[, the (unsealed) 21 entry (e5'^r/2,4)(C^^-' , ?7^^^) of -f^Ar/2,4 satisfies 



(2.1) 



2(e5^/2,4)(f(^\r;W)- 



00 



KAiry{t,ri)dt 

Ai{t)dt){ / Ai{t)dt) 



< o(l)e-'=«e-"'' 



uniformly for ^, 77 G [Lq, +00). This completes the proof of Theorem II . II for /3 — A. 

In the case /3 = 1, we use (|1.42ll and again consider the CD part and the correc- 
tion term separately. The properly scaled 11 and 22 entries of -ftrj^i"'''(C'^'i '7^^') 
converge to the corresponding entries in H1.27|l et seq. with the error estimates 
0(1)6"^^^ and o{l)e~'^^, respectively, uniformly for ^,77 € [Lo,+oo): this follows 
from (|3.8|) for the CD kernel part (giving rise to a smaller error o{l)e~''^e~'^^) and 
from (|4.49|) for the correction term. The properly scaled 12 entry converges to 
the corresponding entry in (|1.27|) et seq. with error o{l)e~'^^e~'^'^ , uniformly for 
^, ?7 € [Lq, +00): this follows from H3.8|l for the CD kernel part and from H4.39|l for 
the correction term. Finally, in view of H1.49|l . H3.56|l and 1(4. 50(1 . the (unsealed) 21 
entry of K^^[''^\^(^\r](^'>) satisfies 



(2.2) 



(e^^,i)(e(^\^(^') 



Ai(s) ds 



K Any {tl])dt 



Ai{t) dt 



Ai(s) ds 



< o(l)e-"™"«''') =0(1) 
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with the uniform estimate o(l) for ^,7] > Lq. In order to obtain the same form for 
the hmit as claimed in Theorem ll.il we note that for aU ^, 77 G M 



/V If f°° 

KAiry{t,v)dt+-{J Ai{t)dt){J Ai{t)dt) 



(2.3) 

KAiry{t,7])dt+-{ Ai{t)dt){ Ai{t)dt). 

^ h h 

Indeed, a direct calculation using the representation H1.25|l for i^^Airy shows that 
the RHS of (|2.3|l is skew symmetric in ^ and r\. In particular, the RHS vanishes for 
^ — 1], as is also evident for the LHS. But the ^ derivatives of both sides are equal 
and hence the identity follows. This finishes the proof of Theorem ll.il 

2.2. Proof of Corollarv ll.3l The following basic fact is well-known (see e.g. |ReSip . 

Let D ^ d/dx denote differentiation and let p{x) be any positive function such that 
e L^{R). Then the operator 

(2.4) A^-^— 
^ ' pD+I 

is Hilbert-Schmidt in L^(M). Indeed, by the Fourier transform, A is unitarily 
equivalent to an operator with square integrable kernel {p~^){k—k')jjj^-^, k, k' G M. 

2.2.1. The case (3 — 2. Let Ai denote the largest eigenvalue of the matrix M in the 
unitary ensemble. It is well-known (sec e.g. ^TW2. ) that for finite N 

Lo 



ProbjAi < cn{i 



In 



= detfl-^^if^(fW,r7W) V 

Since Kn is finite rank, it is indeed trace class. As the trace class determinant is 
continuous under the trace class convergence, we only have to prove that 

(2.5) ANii,v)^^^^§y^KN{e''\v^''^)-KMryitv)^0, aS iV OO, 

in the trace norm in L^([Lo, -l-oo)), in order to prove Corollarv ll.3l for (3 — 2. Let 
xfj,(0 be a C°° function such that xt^iC) = 1 for ^ > Lq and xt^i^) = for 
^ < io — 1- We will show that 

(2.6) iV^oo, 

in the trace norm in L2(M). But then xlo^nXl„ = XLo{xto^Nxta)xLa also 
converges to zero in trace norm in L^(R), where xlo is the characteristic function 
of [Lq, +00), and this clearly proves (|2.5|) . 
Let p(^) = (1 + ^2)1/2 and write 



xt^Nxt = [-p-^] [{D + I)pxlANxt 



*] 

oJ 



The first operator is Hilbert-Schmidt (see (|2.4|) 'l and the second operator is of order 
0(7V~2/3) in Hilbert-Schmidt norm by H3.8|l . with Lq replaced with Lq — 1. This 
proves H2.6|l . 
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2.2.2. The case (3 — A. Let Ai denote the largest eigenvalue of the matrix M in the 
symplectic ensemble. Then in |TW2| the authors prove 



ProbjAi < cn{i 
= Y^det {l 



aAr7V2/3) 



IN 



L2([Lo. + oo)) 



The proof will therefore be complete if we could prove that all the four entries of 
4 (^'■^'' ' ) converge to the corresponding entries of K^^^ ?/) in trace class 
norm in L'^{[Lo,oo)). Again we use p.48|l and prove the trace class convergence 
of the CD part and of the correction term separately. The trace class convergence 

together with the trace class convergence method in Subsection 12.2.11 

To prove the convergence in trace class for the 11 and 22 correction terms, we 
must show that 



of the CD parts of all the four entries of K]^)^'^ follows by using H3.8|) and (|3.5t)|) 



cn 



'f2{e''Y-D21 



1 



D2iCi^ BiiDi2 



^2{t)dt 



Ai(0 / M{t)dt 



(cf. (|4.18|l . H4.22|l ) converges to zero in trace class in ^^([Lo, cxd)). But A^r is an 
operator with finite rank at most n + 1 = 2m = deg V , independent of N . For such 
operators we have the following inequality 



(2.7) 

where 11 • 11 1, 



||Ajv||l < V2^||AAr|jffs 

\hs denote the trace norm, Hilbert-Schmidt norm in L'^{[Lq, oo)). 



respectively. Indeed, |Ajv| = ^^A^Aat is also an operator of rank at most 2m, and 
hence it has at most 2m nonzero eigenvalues, cri > (T2 > • ■ ■ > cr^ > 0, < < 2m. 
Thus 

J \ 1/2 



lA 



trlA 



3 

E 

1=1 



i=l 



< V2m||Ajv|Us. 



But from HTm . \\An\\hs = o{l){J^^ J^^ e-^^e-'"" d^dr]Y^^ = o(l), iV -> c5o, and 
we conclude that ||Ajv||i 0, N oo, as desired. A similar argument using (|4.17|) 
for the 12 entry and (|4.26|) for the 21 entry, completes the proof of CoroUarv 11.31 
for ^ = 4. 

2.2.3. The case (3 — 1. Let Ai denote the largest eigenvalue of the matrix M in the 

3(0 



orthogonal ensemble. Let5(^) = ■\/T+^ and set G'(^) 



. Note 







that 5-1(0 e L\^). Let g^N){t) = ^ I + [^^^{t - a, - dN)? and G(n){0 
9(N){0 

Note that g{N)i£,'"^^) — giO- Recall the definition of det2 

9mi0 ' 
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in the Introduction. A shght modification of the calculations in |TW2I Section 9] 

shows that 

(2.8) 

Prob{A,<c.(l + -^)+rf4 



y det2 ( 



^ _ Cat 



L2([Lo,+oo)) 



In |TW2I Section 9] the authors use the fact that det(l + AB) = det(l + BA) 
for appropriate operators A and B. But one clearly has the freedom to write 
AB = AG'^^jG^N)B, and so we also have det(l + AB) = det(l + AG'^^^G(n)B) = 

det(l + G (jM)B AG~^^^) and this leads to (|2.8|l . We have chosen G(7v) as above in 
such a way as to ensure that 1 + G {^jsi)B AG^^^ has a 2-determinant, but there is 
clearly great freedom in the choice of g[N), and hence of G(jv)- From H2.8|l we see 
that in order to prove H1.33(l it is enough to show [SI| that the diagonal (respectively 
the off-diagonal) entries of ^/j^^/a (G(Af)Jt:^/'^^^G^j^))(gW,7?W) converge to the 

respective entries of {GK'^^^G~^){£,, rf) in trace (respectively Hilbert-Schmidt) norm 
in i2([Lo,oo)). 

We consider first the 11 entry (again the 22 entry can be considered similarly). 
This entry has the form 

where Sn,i is given by the CD part and the correction term as in 11.42|l . The proof 
that 5(0 [-J^Xw(^W,ryW) - KMry{£,.n)\9-\'n) ^ 0, iV ^ oo, in trace norm 
in L^([Lo,oo)) is completely analogous to the /3 = 2 case in Subsection 12. 2. II fnote 
that g and its derivative are polynomially bounded) and the details are left to the 
reader. 

As in the /3 = 4 case above, the fact that the correction term in the 1 1 entry has 
a fixed maximal rank independent of N implies that the trace norm convergence 
follows from the Hilbert-Schmidt convergence. But by 14.40|l . (|4.49|) 



9{0 



<i>i(?(^>)^-Gn- f - r ^,{t)dt\-^,{£,^''Y -G.^-i- r <^2[t)dt 
_ $i(e(^))^ • Gil • e$i(+oo) - $i(e(^^)^ • Gi2 • e$2(+oo) 



1 r 

^ J -oc 



Ai{t) dt 



< 0(1)5(06-^-1(7?) 



which is o(l) in Hilbert-Schmidt norm in L'^{[Lo,(X))). This proves the trace class 
convergence of the 11 (and similarly of the 22) entry. 

Finally, we note from the uniform pointwise bounds in (|1.29() that the error terms 
in the 12 and 21 entries are bounded by o(l)5(0e~'^^e~"'5(ry) and o{l)g^^ {£^)g^^ (rj) , 
respectively, uniformly for ^,7] > Lq. This immediately implies the Hilbert-Schmidt 
convergence of the off-diagonal entries to their appropriate limits. This completes 
the proof of Corollary 

Remark 2.1. With a little more work one can show that in the (3—1 case the 
off-diagonal entries (apart from the term 5-^(0 sgn(^ — 'ri)g~^{r])) in fact converge 
in trace class norm, and not just in Hilbert-Schmidt norm. 
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Remark 2.2. As noted earlier, there is considerable freedom in the choice of the 
auxiliary function g. We see that all we need is that 5, g' are polynomially bounded 
and .g-i e L2(M). 

3. The edge scaling limits of the Christoffel-Darboux {(3 = 2) 
kernel, and of its derivatives and integrals 

3.1. Auxiliary facts from |DKMVZ2j . We now recall some notation from [ibid.]. 
Let dii-^°'\x) denote the equilibrium measure (see e.g. jSaToj ) for OP's correspond- 
ing to the rescaled weight e~^^"(^), Vn = j^V{cnx + (In), where cat, cLn are 
the so-called Mhaskar-Rakhmanov-Saff (MRS) numbers (see |MhSal IRaj ) . For 

V{x) = K2mx'^™ + K2m-ix'^™~^ H as in H1.19|l . we have [ibid., Thm. 2.1] to any 

order g as 00 
(3.1) 

( J }^"')- ] ivi/(2™) +y C(,)A--'"/(2") +0(A^-(''+i)/(2'")) 



cn 
and 



(3.2) dN - -^H!^ + y ^ Ar-J7(2™) + (9(^-(9+i)/(2m))_ 



As Af — !■ 00, the equilibribum measure is absolutely continuous with respect to 



Lebesgue measure, dfj,'-^'^^ (x) — (x) dx, and is supported on the (single) interval 



[-1,1], 

(3.3) i>^P\x) ^ ^n{x) = ^|1 - x^\"\^_^^,^{x) hr,{x) 

(see [ibid., (2.4)]) where h^ix) is a real polynomial of degree 2m — 2 satisfying 
[ibid., Prop. 5.3] 

(3.4) /lAr(x) > /imin > 0, X £ R, N>Ni{V). 

Set [ibid., (5.33)] 

g{z) = gN{z)^ J ^og{z ~ x) dn'-^'^\x) = J log(z-x)^|l-a;^|^/^/iAr(x)dx, 
z e C \ (-00, -1], and for z e (-1, 1) [ibid., (5.34)] 

(3.5) EnIz) = g+{z) - g^{z) = i \1 - x'^l^^'^hNix) dx. 



We also use the same symbol for the analytic continuation of ^jy to C\ ((—00, — 1] U 
[l,+oo)). 

Notational remark: Here we denote by what was denoted by £,n in [ibid.]. 

For a fixed 6 > sufficiently small (cf. |DGI Rem. 4.3]), let R denote the matrix 
function defined in |DKMV!Z2l (7.47)]. The function R is analytic in the comple- 
ment of the contour S/j as in [ibid., Fig. 7.6] and is continuous up to the boundary. 
Furthermore by [ibid., Thm. 7.10], it has an asymptotic expansion 

00 

(3.6) i?(z) -/ + A-i^rfe(z)A^-'=/(2™) 

fc=0 
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where {rk{z)} are bounded functions that arc analytic in the complement of {\z — 
1| = (5} U {jz + 1| = 6}. The expansion (|3.6|l is uniform for z G C \ S^. Moreover, 
by the proof of [ibid., Thm 7.10] and Cauchy's theorem, it follows that H3.6|) can 
be differentiated term by term, 

(3.7) ^R{z) - f2 ^r,(z)JV-'=/(2™), j = 1, 2, • • • , 

k=Q 

where again the expansion is uniform for z e C\Efl. Also, each -j^fkiz) is bounded 
(and analytic) in the complement of {\z — 1| = (5} U {\z + 1| = 5}- 

3.2. Estimates on the CD kernel and its derivatives. We will only consider 
the end point 1 (the end point —1 can be considered similarly). Let io G K be 
fixed. Recall the notation in H1.25|l . (|1.23|) and (|1.5|l . Our goal in this Subsection 
is to prove that for j. A: = 0, 1, and some C — C{Lq) > 0, c = c{Lq) > 0, one has 
uniformly for ^,1] Cz [Lq, +oo) 



(3.8) 



Cn 



rif^(e^^\r7^^')-ifAiry(e,ry) 



Note that [AbStj 

I Ai(e)| < Cil + |e|)-^/^ I Ai'(ai < Cil + 1^1)1/4, e e K, 
|d«Ai(OM'| <Cie-«<C2, Ce[Lo,+oo), <7 = 0,1,2,--- . 

3.2.1. Auxiliary notation. Set (see |DKMVZ2I (2.15)] and also DG, (4.10)]) 

(3.10) fN{x) = aw7V2/3 {x - l)fN{x) 

which satisfies the following (see (the proof of) |DKMVZ2l Proposition 7.3]) 

(1) /jv(a;) is real analytic on (1 — 2(5, 1 + 2(5), and to any order (7 = 0, 1, 2, • • • 

fN{x) ^ ^iV-J/(2™) /(^.)(x) + 0(Ar-(9+i)/(2™)) 
i=o 

uniformly for x in the interval. Moreover, the functions /(j) (x) are also real 
analytic on 1 — 2(5 < x < 1 + 2(5 

(2) to any order q — 1,2, - ■ ■ 

aN = (/i?,(l)/2)'/' = 2to2/3 + J2 iV-j/(2m) ^^^^ _^ Q^^_(,+ l)/(2™)^ 

i=i 

(3) f'j^ix) = -ajviV^/^VFAr(x), where Wn{x) = /Ar(x) + (x - l)/;/(x) also has 
an expansion uniform in x to any order q = 0,1,2, ■■ ■ as above 

q 

3=0 

The terms (x) are real analytic on 1 — 25 < x < 1 + 25 

(4) maxfc=o,i:2maxi_25<2:<i+25 \d'^ fN{x) /dx''\ < M < 00 for N > N2(V) 

(5) /^v(l) = 1 = Wn{1) and mini_25<.<i+25 /iv(x) > i for > Ar2(T^). Also 
/(o)(l) = l = VK(o)(l). 
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Denote 
(3.11) 



= e/(«JviV2/3), r,N = 7i/{aNN^h 

In = [Lo, onN^^^S], I In = [onN^^^S, +oo). 



Thus, recalling ifT"^ . = CN{l+^N)+dN and similarly 77^^) = CN{l+VN)+dN- 
As above, let 5 > be fixed and sufficiently small. Consider ffi'st ^NtVn in a 
neighborhood of 0. Set 



(3.12) 

and also 
(3.13) 



9N{0 = ^fN{l + ^N) 
Fn{1+(n) = + (/Ar(l + eA'))'^' 

Fn{1+Cn) = N^/''a]I^FN{l + ^N) 
AoiO = N'^'^a'j^^FNil + (n) ■ AiigNiO) 



= /v-1/6^71/4 



FNii + i 



Note that in view of H3.10|l fl')f5'l and the formula 

(3.14) 9'N{0^fN{l+^N)+^Nf'N{l + ^N) 

there exist C2 > ci > such that 
(3.15) 
and 

(3.16) ci < g'NiO < C2, WnIOI < C2N-'/\ 



Cl < < C2, 



N 



Fn{z)\ 



< C{k) 



Similarly one has uniformly for £, & In 
(3.17) ci<Fn{1 + ^n) <C2, 

for some C{k), fc = 1, 2, • • • . 

3.2.2. Estimates for (^, rj) In x In ■ With the above notation the following holds. 
Proposition 3.1. For (C,??) £ In x In 

Cat ^, ,_.,An tAn-. s 1 



(3.18) 

where 



if^(^W,,7W)=i^W^,r;) + 



n,i 



1 -i 

-1 -i 



(3.19) 



{R^y{i + ^N+t{m-iN)) dt 

-1 



1 -i 

-1 -i 



-MO 



18 



DEIFT AND GIOEV 



and 
(3.20) 



lo ^-iv(l + VN + t{^n - Vn)) dr 



and 



i^Ar(l + w) 



(3.21) 



Qi,3i^,r^) = ENitv) / Ai(z + 5jv(0)Ai(z + ffw(r?))dz 



+ + <^{VN + t{^n ~ w))) do 

Jo 



dr 



(3.22) 



i,4(C, V) = ^LNiO / C^w(^, ^) Ai(z + .gw(77)) 



+ tj^Ln{7i) / Ai(z + OC/w(77,-2)c;-z 

Jo 



Proof. First, some algebra: let Y solve the Fokas-Its-Kitaev Riemann-Hilbert 
problem for the polynomials orthogonal with respect to the weight e'^^'^^dx (see 
IDKMVZ2I Thm. 3.1]). Then as in (DKMVZTI (6.3)] we find for any x,y eR 



(3.23) 



KN{x,y) = e-(^(")+^(«»/2 Yuiy)Y2iix)^Ynix)Y2iiy) 

^ __^_(y(.)+v(,))/2 ( 1 0)-Yjiy)-Y^^{x)-{0 1 )' 

2'Ki{x — y) 



Here and below +/— refer to the boundary values taken from above/below the 
real axis. (The choice F+ is made only for definitcness. Formula H3.23|l clearly 
remains true if Y-^ is replaced with F_.) Consider first z = 1 + ^at € (1 — (5, 1] for 
C G (-<5aAriV2/3,0]. By |bRMV^2l (4.2), (4.6), (4.22)] we have for 5, the solution 
of the Riemann-Hilbert problem [ibid., (4.24)-(4.26)], (cf. [ibid., (7.46), (7.47)]) 



S+{z) =c 



■N 



y+(cArZ + d 



(3.24) 



X e 



-A(s+(2:)-l)<T3 
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where (73 = and the constant I = In is given by [ibid., (5.35)]. Solving 

Vo -1; 

for F+ and substituting in (|3.23(l we find for ^,r] £ {—6aNN^^^,0] 
CN (N) (N) e--T(^"(i+«")+^"(i+')")) 



(3.25) 



X ^ gA'(s+(l+7)«)-^) gV(g_(l+,,„)-i) ^ .sl{l + r]N) 



=JV(9-(l+?iv)-^) 



.JV(s+(i+5„)-i) 



Now note that for z G (1 - 5,1], by [ibid., (7.46), (7.47)], S{z) = R{z)Pn{z). By 
[ibid., (7.24), (7.9), (7.23), (7.4)], 

1 -1 \ / Fn{z) 



Pn,+ {z) = y/ne 



ITT /6 



1/Fn{z) 

( 1 M 

X ^/+(/Ar(z))e-"'^-^/6 ( I ^N^M{z)a^/2 

where 

/ Ai(/^(z)) Ai{io^jN{z)) \ 
AI{fN{z))^\ \, c. = e2"/3. 

V Ai'(/^(z)) iv^Ai'{u^fN{z)) j 

For z G (-1, 1), in view of [ibid., (5.38)] 

-VN{z)+g+{z)+g^{z)-l^O 

and we find from (|3.25|l 
(3.26) 



anN-^'^ 27ri(e - ?/) 

^^w(l + w) \ I ^ 

l/F^(f + r,Ar) / V -1 -I 



-1 



X i?;^(i + ^A,) 



1 -^ \ / l/Fjv(l + ew) 



-1 / V Fn{1 + ^n) 

X ^/+^(/Ar(l + Cat)) j , e, ^ e (-<5aJviV'/^ 0]. 

Similar calculations for z G [1,1 + 5) lead to the same formula for all other cases 
C < 0, r; > 0, etc., Id, |?7| < fojviV^/^- 
Now writing 

(3.27) R^{\ + m) = ^^(1 + Cjv) + {m - f a) / (i?^)'(l + ^A + i(w - f a)) dt 
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and taking into account dctyl/+(/Ar(l + £,n)) = -l/{2T:ie"/^) (use [ibid., (8.38)]) 
we obtain from (|3.26|l that 

A^^(eW,^W) = ^_Q,,,(e,,) 



(3.28) 



t - ?? L ^iv(l + VN) 

where Qi.i is as in H3.19(l . Now 

^3 29) Fn{1 + Vn) Fn{1+vn) 

= 1 + {^N- Vn) — ^— ^ 

Fn {1 + m) 

and hence using (|1.25(l we rewrite (|3.28(l as 

ir (c(N) (N)\ 9n{£.) - 9n{ii) ( (t\ I \\ 

^33Q^ ^^^^^(e^ K^..y{9N{^UN{v)) 

where Qi,2 is as in Next we write ^"^^^"^"^''^ = 5^(77 + t(^ - ry)) dr, and 

use (|3.14ll and 

/w(l+'77V+T(CAf-'7Af)) = /w(l) + (?77V+T(?Ar-?7Ar)) / (l + CT(j?Ar+r(^Ar-77Ar))) dcr 

Jo 

to conclude that ^"^^^"^"^''^ = l + £:jv(^, r;) from (E^U- Hence recaUing ((T^ we 
obtain from (|3.30() 

(^•31) ^^^-(^'"^'?^"^)-^A.y(5A^(0,5A'(^))+--^ 

where Qi^a is as in l|3.21|l . Finally again using (|1.25|) we find 
(3.32) 

KAiryigN{0,gNiv)) = / Ai(z + 5Ar(f))Ai(z + gjv(»7))c?2: 



Ai(z + ^) Ai(z + 77) dz 

+ / Ai(z + 0[Ai(2; + 5Ar(77)) - Ai(z + 77)(i2; 
Jo 

/•oc 

+ / [ Ai(2: + 5Ar(0) - Ai(z + 0] Ai(z + gNiv)) dz. 



The first integral equals KAhyi^, v)- To evaluate the third integral we recall H3.12|l . 
()3.22|l and note that 

9N{0-^^^[fN{l + (N)- fN{l) 
(3.33) ^2 rl ^ ^2 
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which impHes 

Ai{z + QNiO) ~ Ai(z + e) = ^J^y, ■ ^ /^(l + <Jv) da 

(3-34) X f\i'{z + ^ + T{gN{0-0)dT 

Ja 

The second integral in (|3.32|) is treated analogously. We conclude from H3.31|l . 
(I3.32|l that H3.18|l holds where Qia is as in (|3.22l) . The proof of Proposition 13. II is 
complete. □ 

Now we prove the estimate H3.8|l for ^,77 £ Ijy. Note that by H3.15|l it follows 
that gNiO^ QNiv) bounded below by some constant A/q, and hence in the region 
{^jf]) In ^ In, both variables are bounded below by the constant Lq. Using in 
addition H3.16|l we conclude that we can always use the exponenial bounds on Ai 
and its derivatives in H3.9I) . and hence for any m £ N and fc = 0, 1, 2, as ^ 00 

(3.35) |^™(^)^-Ai(g^(^))| <C(m)e-^(™)«, ^eIn- 

Consider Qi, 1(^,77) first. Recall from H3.7|) that, in particular, j^R{z) — 0{N^^), 
j = 1, 2, 3, uniformly for |z - 1| < 6. It follows then by using ifTT^ . ifTTTjl . 

(IX^ that for j, fc = 0, 1 

(3.36) |9|5^'Qi,i(C,77)| < const • iV-^'/^g-cCg-c^ 

uniformly for ^, i] £ In- In the same way we find that for j, fc = 0, 1 and ^ = 2, 3, 4 

(3.37) |9|9^Qi,i(e,ry)| < const • iY-^/Sg-^^g-"' 
uniformly for G In- In estimating Qi,4, we use the estimate 

\9NiO-^\<Cm\, \i\<SaNN^/\ 

which follows from (|3.33l) . together with the uniform boundedness of Ln{£,) (see 
(|3.1U|l f4')'): for S sufficiently small this implies that 

(3.38) \UNii, z)\ < Ce-'^e-'^, ^eIn, z > 0, 

with similar estimates for the ^- (and z-) derivatives. This proves (|3.8|l for (^, rj) e 
In X In- 

3.2.3. Estimates for (^,77) £ IIn x IIn- Recall from [ibid., (4.30), (4.31), (6.16)] 

a(z) + a{z)^^ i{a{z)^^ — a(z)) 



(3.39) S^°°\z) = N{z) 



i{a{z) ^ a{z) ^) a{z) + a{z)' 



where a{z) = (fqrf)^^* 1 as z — > oo- 
Proposition 3.2. For j, A: = 0, 1 and some C,c > 



(3.40) 



uniformly for £,,ri £ IIn- 
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Proof. Note first of all that (|a.23|l still holds. For z = 1 + ^ [l + S, +00) we now 
have in place of H3.24|l 

(3.41) S+iz) = c-/''' e-^"^ Y+{cNZ + djv)e-^(s+(")-^)"3 

where again the constant ^ = Zjv is given by [ibid., (5.35)]. Solving for Y+ and 
substituting in H3.23|l we find for ^, 77 G IIn 
(3.42) 

./^j^(^W^^(JV)-) ^ _ g-f(y„(l+e„)-23+(l+5„)+0g-f (VV(l+,,„)-23+(l+,,„)+i) 



Cn 



( 1 ) •5j(l + 77jv)5+^(l+eiv)- ( 1 )' 



2m{^-vi) 

In view of [ibid., (5.38)] 

-Vn{1 + £.n) + 2.g+(l +iN)-l^ SAr,+(l + ew), e e //w. 

Now by [ibid., (2.14), (5.34)] for some Ci((5), C2((5) > and c > for N large 
enough 

< - / V^^-l^min - / c/lmin 

< -Ci - C2(C7V - <^), ^ e (foAriV2/3, +cx3). 

By [ibid., (7.46), (7.47)] for z > 1 + 5, S{z) = i?(z)S'(°°)(z). Using (TTTTIl . which is 

still valid for ^, 77 € //at, we obtain 

(3.44) 

5j(l + w)5+^(l + Ca) = 4°°^'^(1 + w)5l^)-^(l + Ca) 

+ im - CN)si°^^^ii + w) Q (i?^ )'(i + Ca + i(r/A - Ca)) d?j s^_^^-^ii + ^n). 

Substituting 

+ ivN (^i°°^^)'(i + ^N + t{m - ^n)) dt 

in the first term in the RHS of and noting that (l 0)-/-(0 1)^ = 0, 

we obtain an expression for ( 1 ) • 5^°°^^(1 + i]n)S^+°^ "^(1 + ^jv) ■ ( 1 ) 
which is proportional to (^ — r/). The exponential bounds H3.40|l then follow from 
(IS^ and the properties of S^°°^ and R (see and respectively). □ 

Now we prove (|3.8|l for 11^ by showing that both of the two terms on 

the LHS of H3.8(l satisfy the exponential bound. More precisely, let ^ G I In- Then 
either ^jv > 2^ or g [5, 2(5]. In the former case 

(3.45) e-=^«"-*) = g-cA((u/2)-5)g-(cA/(2.„AV3))^ ^ TV 00, 

since — > (2m)^/'^ ^ as ^ 00. In the latter case 
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and hence 

(3.46) e-^^e-^^«"-*) < e-(-/2)^e-(^/^)^ < e-(^/2)^e-«, e [S, 26]. 

Combining H3.45|l and (|3.46|l we conclude that Proposition 13 . 21 implies 

Now we consider KAiry{£,, v) for ^, 77 G IIn- It follows from jAbStj that 

(3.48) I Ai(a;)|, | Ai'(a;)| < C{Lo)e-<^o)M'''\ x > Lq. 
Using the integral representation (jl.25|l we estimate for r/ G IIn 

(3.49) |ifAiry(e, V)\<C e--(-+«)^^%-c(^+.)-^^^ dz. 

Jo 

Let ^ e IIn- Then ^ > 1 for large N. It is elementary to verify that 

(3.50) {z + £,f/^>z^^^ + £,^/^, z>0, ^>1. 
Next, ^3/2 > ((5aAr)^/2^ > -> 00 and hence 

(3.51) ^^/'^ ~^ = ^^^\l~C^^^)>c"N, N-^00. 
Inserting (|3.5Q|I . I|3.51|l and their analogs for 77 in (|3.49(l we find 

(3.52) \KAnyi^,v)\ < Ce-^^e"^«e~-', C,7j e IIn- 

A similar argument using 13.48|l also shows that the derivatives of if Airy satisfy 
the same bound. Combining 13.47|l and H3.52|l completes the proof of (|3.8|l for 
IIn- 

3.2.4. The "mixed" neighborhoods of the end point I: (^,77) € {In x I In) U {IIn x 
In)- Let us consider the case (^,r?) G In x I In (the other case is treated analo- 
gously). For i^Airy, we use the bound in H3.9|l for ^, 

|Ai(z + OUAi'(z + OI <C(Lo)e-^e-«, z > 0, ^ > ^0, 

together with the bound (|3.48|) for 77. Inserting these bounds in p.25|l we obtain 
for j, A: = 0, 1 

(3.53) \dld^KMry{tv)\<Ce~^''e~'^e'^\ (^, 77) e x //^ 

as before. For ifAr(^(^\ 77^^^), there are two cases: |^7v — wl ^ ^"^^ — 
r]N\ > S/2. In the first case we can treat both points as lying in a /at x In region 
corresponding to a larger (fixed) value of S (more precisely, set 6 36/2) and hence 
(|3.8|) follows by the arguments in Subsection 13.2.21 

It remains to consider the case {£,, 77) e /at x IIn, \£,n — i1n\> S/2. For such ^, 77, 
we have 

(3.54) 1^ - 771-1 < N-^/^a^'^26-\ 



(3.47) 



Cn 



« naA^2/3 



ifA(e^^\r7W) 
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The computations that led to l|3.26(l and (|3.42|l now imply for E In jf] ^ 

Cm /.n g-7ri/6 -ArHjv(l+U)/2 



(3.55) 



X ( 1 ) • 5(°°)^(1 + m)Rl{i + Vn)R+^{1 + ^n) 
1 



-I ~i 

Xyl/-^(/^(l+eAr)) 



l/i^w(l + eiv) 

Fn^l+iN) 



1 



Kn{s, 77^^)) ds-{- \ KAny{t, v) dt) 

(N) ^ ^ 

ifjv(s, 77^"^^) - ( - / -^Airy(i, r;) 

(«) ^ Je ^ 



and using the preceding estimates we find for j, fc = 0, 1 

uniformly for (^,77) g /at x //at, I^at - wI > <5/2- 

There is a similar estimate for (^,77) € //jv x which, together with H3.53|l . 
then proves H3.8() for (^,77) e (/at x 11^) U (//at x /at). 

3.3. Estimates on integrals of the CD kernel. For ^,776 [L07+00), making a 
change of variables s — cn{1 +1^) + d^ , and using (|3.8|) with j = fc = 0, we readily 
find 
(3.56) 

poo ^ poo ^ 

uniformly for ^, 77 G [Lqi +00). 

4. The contribution of the correction term for /? = 1 and 4 

4.1. Auxiliary facts concerning integrals of the orthogonal functions 0^. 

It was shown in jPCxI (4.14)] that for a fixed G N the following holds as ^ 00 
(see (TTTni 'i 

(4.1) / 0^+,(y)d7/ = c]^/^Ar-i/2(2m)-V2(i + (_i)A^+i + O(Ar-V(2™))) 



where 277z — deg V^. Introduce the following column vectors of size 2m — 1 
a=(l,0,l,0,--- ,1)^, b=(0,l,0,l,--- ,0)^ 
e = a + b = (1,1,1, ••• ,1)^. 

By (|1.47|l and (|4.1|l as (even) N ^ 00 



(4.3) 



e$i(+^) = - / $1(77) dT/ = c]^/' (2777)^1/2 (b + 0(1)) 



' —00 

/' + 00 



e$2(+oo) - - y $2(y) = c]^^' (2777)-i/2 (a + 0(1)). 

We need also the following result. Recall the notation H1.23|l . H3.11|l . 
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Proposition 4.1. For any fixed j G N there exist C, c > such that the following 

holds as N —I- oo, 



(4.4) 



<t>N+At^ ^) - ^ 1/2 



This estimate implies that for a fixed j (z Zi there exist C, c > such that 
(4.5) 



i-oo ^1/2 2l/4 (-oo 



Ai(i) 



i V 

Proof. Assume first that j = 0. It was shown in | DKMVZ2I Thm. 2.2] that (in our 
notation) 
(4.6) 



(t)N{t 



N 



a]^^N^/''FN{l + tAr) Ai{gN(t)) (l + 0(7V-i)) 

^-i/4^-i/6(^^(^ + iAr))-i Ai'(5^(i)) (1 + O(iV-i)) 



where the error terms are uniform for t In. Using (|3.35() we immediately estimate 
the second term above by Cc^^"^ N~^l'^e^^^ uniformly for t G 1^. The part of the 
first term that corresponds to 0(iV~^) is estimated similarly by Ccjy^^^A^^^/^e"'^*, 
t £ In. To estimate the leading part of the first term we write 

Fn{1 + tN) Ai{9N{t)) = i^Ar(l) Ai(t) 

+ Fjv(l)[ Ai(.9Ar(t)) - Ai{t)] + AiigNit))[FN{l + tN) - Fw(l)] • 
By (jSIIKS) and (Tm . Fn{1) = 2^4. By formula iiOHl and XTM 

I Ai{gNit)) - Ai{t)\ < CN-^^H^e~'' < CiV-^/Sg-c't^ ^ ^ 
Also using H3.35|l and H3.17|) we obtain 

I Ai{gN{t))\ ■ \Fn{1 + tN) - Fn{1)\ < CiV-2/3|i|e-ct < c' N-^/^e-'"'* , t G In- 
Combining the above estimates we find that 
\FN{l + tN)Ai{gN{t))-2'/^Ai{t) 



e I 



which completes the proof of H4.4() for j = and t £ In- We now consider 1)4. 4|l for 
j = and ^ e IIn = [SonN'^^^, oo). For such t, by lpO?|l 



|Ai(i)| <Ce-*<Ce-^^'''e-*/= 



and hence |^^^^^^^^^^ Ai(t)| < Cc^^^^iV-i/'^e-^*. Also from (4.8)], we find 

|</'Ar(t*-''^-')| < Cc^^^^e~'^^e~^*. These two estimates for t g IIn, together with the 
previous estimate for t & In, yield (|4.4|) in the case j = for alH G [Lq, oo). 

Now fix any j £ Z and write 
(4.7) 



4.7) 
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where 



(4.8) 



tN./ 



t ■ 



CN aN+j {N + / CAT 



CN+j O-N 



CN+j \dN+j 



iV2/3 



\CN+j 



= (l + 0(7V-i))t + 0(iV-i/3) 

by (EH), lE!2l, (|XTn|) f2). In particular, as TV ^ oo, ijvj > (1 - 5SgnLo)-^o- Now 
the RHS of (|4.7I) can be written as (/)jv' ((ijv,i)^^'^) where N' ^ N + j. Applying 
the estimate (|4.4ll just derived for j = 0, with Lq replaced by (1 — ^ sgnio)-^o, we 
obtain 
(4.9) 

c^NiitN.r'^) - ""'^•^'^,)/^''''''' Ai(t.,-) < Cc-/;(iV + ,)-/^e--- 

for all t > Lq. Using (|4.8|) . and also H3.1|l . (|3.10() ('2'). together with the elementary 
estimate 

I Ai{tN,j) - Ai{t)\ < CiV-i/^e""'* 



□ 



(use (j^J), we obtain (|4.4(l from (|4.9(l for any fixed j e Z. 
Finally 14.5|l follows readily by integrating H4.4I) . 

Recall the notation (|4.2|l . Proposition 14.11 implies that for j = 1,2 one has 
uniformly for t,^,r] > Lq 



$,(tW) 



1/4^1/621/4 

^ ■ Ai(t) e 



.1/2 

-TV 



(4.10) 



/ *j(s)'^s 



Ai{t)dt]e 



,(N) 



4.2. The case /3 = 4. 



/ 4' 21/4 



Ai(t)di 



4.2.1. T/ie contribution of the correction term to the 12 entry of K^a- Since {SD){x, y) 
—dyS{x,y), the correction term in p.43|l has the form 

(4.11) -$2(x)^ • D21 ■ $1(2/) - <i>2{xf ■ D2iC^^^BnDi2 ■ Mv)- 

Set X = y = ryf^). Recall n 2m- 1, 2m = degF. Note that by jDHI (2.13)] 



(4.12) 



mf^2m 2m-l 



22" 



-1 "-TV 



/ (0) 
1 



V 



{(n-l)/2}\ 



1 / 



0(1) 



and D21 = 0(c^""^) = 0(Ai-i/(2m)) as A^ ^ 00 (see (|XTJ). Also since Cn = 
/ — B12D21, we see from [ibid., (2.19)] that C^^Bn is skew symmetric, being the 
lower right corner of the skew symmetric matrix -D^^. (Note that Bn is the lower 
right n X n corner of ejv.) Hence D21C11B11D12 in H4.11|l is also skew, and using 
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[ibid., (2.13)] and the fact that C^Y is bounded as N oo [ibid., Thm. 2.4, 2.6], we 
see that D2iC^iBiiDi2 = 0(iVi-i/(2m)) as oo. Recall that the 12 entry in 

Kna{£,^'^\v^^^) J^^s an overall scahng factor (^^-^273)^- Substituting the leading 
term in the representation of <I>j in H4.1U|I into the first term in H4.11|l . and using 
(|4.12|l . we obtain 

(4-13) ^ (S„ + o(l)) Am Ai(,) 



where o(l) is independent of ^, 77 and S„ denotes the sum of all elements of the first 
(binomial) matrix on the RHS in H4.12(l . Using the formula preceding [ibid., (6.7)] 
one finds 

1 TO(2m)! 
" " 2 (to!)2 ■ 

Recall that by (|XTUI) f2). [ibid., (2.14)] and 

(4.14) r V ^ ^ 22™-i (2m)! ^ 

|Ai(0| <Ce-«, ^>Lo. 

Inserting these estimates, ()4.13|l becomes 

-iAi(0 Ai{v) + o{l)e-^e-^ 

uniformly for f , 77 > io and o(l) is independent of f , 77. The error that was made by 
substituting only the leading term in (|4.1U|) in the first term in (|4.13l) , is estimated 
as follows: 
(4.15) 

0{c%N-^/^) 0(c^™-i) {0{N^/'c-'/')c-'/'N'^/^{\ ^(016-^" + I Ai(r;)|e-^«) 

< 0(iV-i/3)e-'=«e-"'' 

uniformly for ^, 77 > Lq, and independent of the degree 27n of V. 

Next we substitute the leading terms in the representation of $2 h^ H4.10|l in the 
second term in (|4.11|l . By the skew symmetry of D21C11B11D12 noted above, the 
result is precisely zero. The error that is made by such a substitution is estimated 
in exactly the same way as in (|4.15f) and is also of order 



(4.16) 0{N 



-l/3V-cCg-c^ 



uniformly for ^,77 > Lq. We conclude that the contribution of the correction term 
to the 12 entry is given by 

(4.17) -iAi(e)Ai(77)+o(l)e-^«e-=^ 
uniformly for ^,77 > Lq. 
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4.2.2. The contribution of the correction term to the 11 and 22 entries of Kna- 
We consider the 11 entry of Kn,4 (the 22 entry is analyzed in the same way). The 
correstion term in H1.48|l has the form 



(4.18) 



Mxf-D2i-(^-J^ ^i{t)dt 

+ <^>2{xf ■ D2iC^i^BnDi2 • ( - _^ *2(t) dt^ . 

We set X = y = r?^^) in The 11 (and 22) entry in KNA{£,^^\'n^'^^) 

has an overall scaling factor ^^^2/3 ■ Hence, substituting the leading terms in the 
representation of $2, / 'i'l in H4.10|) into the first term in H4.18|l and using H4.12|l . 
we obtain 
(4.19) 

_-^|-.-.fv^jA_,.„,„a))AiK,(-/ Ai(„..) 

where o(l) is independent of ^,7/. Computing the factor and using H4.14|l as above 
we see that (j4.19|l becomes 

--Ai(0 / A{{t)dt + o{l)e-^e-'^ 
2 Jn 

uniformly for ^,7/ > Lq and o(l) is independent of ^,77. The error that was made 
by substituting only the leading terms for <i>2, / $1 in (|4.1U|) into the first term in 
(I4.18|l . is estimated as follows: 
(4.20) 

//Vi/6 1/2 1 1/2 1/2 1 N 

= 0(iVi/3)(iV-2/3 + 7V-2/3 + iv-ije-^^e-^" 
= 0(iV-i/^)e-"«e-"'' 

uniformly for ^,r] > Lq, and again independent of the degree 2m of V. 

Next we substitute the leading terms in the representation of $2, / $2 into (|4.10() 
in the second term in H4.18|l . Again by skew symmetry, the result is precisely zero. 
The error that is made by such a substitution is estimated in exactly the same way 
as in (|4.20|) and also has order 



(4.21) 0{N' 



uniformly for ^,77 > Lq. We conclude that the contribution of the correction term 
to the 11 entry is given by 



(4.22) -^Ai(0 / Aiit) dt + o{l)t 

2 Jri 

uniformly for ^,77 > Lq (for the 22 entry ^ and 77 should be interchanged). 
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4.2.3. The contribution of the correction term to the 21 entry of Km, a- By H1.48|l . 
I)1.5()|l the correction term in the 21 entry of Kn,^ is given by 

/•OO /'OO 

/ <^l{s)ds -021- / ^l{t)dt 

(4.23) 

+ / ^^{s)ds-D2iC^i^BnDi2- <i>2{t)dt. 

•J X J y 

Again we replace x = ^'^^ y = ?7^^'. Recall that the 21 entry in i^Af,4(C'-^\ 77*^^^) 
has no overall scaling factor. Substituting the leading terms in the representation 
of J $j, j = 1, 2, in H4.10|) into the first term in H4.23|l in the same way as before, 
we obtain 

OI/2 /"OO /"OO 

+ / Ai(.)dW m)dt 
^ oij!] N J(, J-n 



ki{s)ds I Ai{t)dt + o{l)e~^e 



uniformly for ^,77 > Lq and o(l) is independent of ^,77. The error just made is 
estimated as follows: 

/ 1/2 1/2 \ 

^/ 2m-lx / o '^N CN \ eg -cr, 

^(Cjv j Ar5/6 + ^5/3 ^ 

= 0(7V-l/3)g-cCg-cr, 

uniformly for ^,?7 > Lqi here all order factors are independent of ^,77. (Here we 
have used | Ai(t) dt\ < Ce~^ uniformly for ^ > Lq.) 

Finally, we substitute the leading terms in the representation of J ^j, j — 1,2, 
in (|4.10() into the second term in (|4.23() . By the skew symmetry the result is again 
precislely zero. The error that is made by such a substitution is estimated in exactly 
the same way as in (I4.24f) and is also of order 

(4.25) 0{N-^/^)e-''^e-''^ 

uniformly for ^,77 > Lq. We conclude that the contribution of the correction term 
to the 21 entry is given by 

-1 />oo noo 

(4.26) -/ Ai{s)ds I Ai(i)dt + o(l)e^"«e~"' 
uniformly for ^, 77 > Lq 



4.3. The case (3=1. As we will see, this case is more involved than the case /? = 4. 
Consider the 2n x 27i (71 = 2to - 1, 277i = degF) matrix {AC{l2n - BAC)~^Y 
in the (3=1 correction term in p.42|l as a two by two block matrix with blocks 
of size n X n. Denote the upper left and the upper right blocks by Gn and G12, 
respectively. With this notation the correction term has the form 

(4.27) -<l>iixf ■ Gn • e$i(y) - $i(x)^ • G12 • 6*2(2;). 



As in |DG| let R = i?„ denote the n x n matrix with all entries zero apart from 
ones on the anti-diagonal (thus Rij ^lifj^n — i + 1, 1 < i < n, and = 
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otherwise). Note that — In. Define 

(4.28) Gil = -i?Z?2iCiVSiiZ?i2i?. 

Note from Subsection lQU that 0210^;^^ B11D12 is skew and of order 0(ivi-i/(2™)) 
as N ^ 00. Hence Gn is also skew and has the same order as ^ 00. We need 
the following result. 

Proposition 4.2. As (even) N ^ 00 we have Gii,Gi2 = 0{N^^^/'^^"^^), more 
precisely 

(4.29) Gil =Gii +o(A^i-i/(2'")), A^^oo, 
and also 

(4.30) G12 = D12 +o(7Vi-i/(2™))^ iV->c3o. 
Proof. It was shown in DG, Theorem 2.3] that, as 00, 

{BA)22=-R{BA)nR + o{l) 





(4.31) 

^ ' BAC = 



(BA)21+0(1) (BA)22+0(1) 



Denote 
(4.32) 

T = In~iBAC)22 ^ In-iBA)22+o{l) = In+R{BA)nR + o{l) ^ RCnR + o{l). 

It was shown in |DGI Theorem 2.6] that, as — > 00, T approaches a constant 
nondegenerate matrix. Thus 



ihn ~ BAC)-^ 



In 

T-\{BA)2i + o{l)) T-i 
and simple algebra using (|1.39() . H1.40|l now shows that in the product AC{l2n 

we have by H4.31|l 




G^^i = A12 [{BA)2i + {BA)22T-\{BA)2i + o(l))] 
GI2 = A21 [In + (S^)ii + {BA)i2T-\{BA)2i + o(l))] . 
Using H4.32(l . this implies 

^-l+l/(2m)^T ^ iV-l + l/(2™)^^2[j^ + {BA)22T-^]{BA)21+0{\) 

(4.34) = A^-i+i/(2™Ui2 [T + {BA)22\T-^{BA)2i + o(l) 

= N-^+^/^'^"''^Ax2T-^{BA)2i + 0(1). 

Now from 

B12A21 i3iiAi2 

(4.35) BA 

B22A21 B21A12 
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and H4.a4;i . (|4.32(l wc obtain 

^-l + l/(2m)(^T ^ Ar-l+l/(2™)^^2(i?i? + R{BA)iiR)-^B22A21 + o(l) 

= ^-i+i/(2'")^i2i?(/„ + {BA)nr'RB22A2i + o(l) 
= N-^+^/^^"'^A,2RC-,^RB22RRA2i + o(l). 
Using the asymptotic relations 

^-l + l/(2m)^^^^^ = iV-l + l/(2™)yl2i + o(l) 

from 'DGI Subsec. 5.2] we see that 

^-l + l/(2m)^gT ^ = iV-l + l/(2")(i?Ai2i?)Cri'(i?B22i?)(i?A2li?0 + o(l) 

= -7V-i+i/(2'")A2iCiVSnAi2 +o(l) 
= iV-i+i/(2™)^2iCiVSni5i2 + o(l). 
As noted above, the matrix -021(^11^511^12 is skew symmetric and hence 

j^-i+i/(2m)fj^_^ = -N-^+^^^^"''>RD2iC^i^BnDi2R + o{l) 



which proves 

Now let us prove H4.30|l . From H4.33|l we derive 

^-i+i/(2m)g,T ^ N'W/i^n.)^^^ ^ (5^)^^ ^ {B A) uT'^B A) 21] + 0(1) 
and hence, because A21 = A12 = -D12, we note that we just have to prove 

^-'+'/^'"'^2i[(Bv4)ii + {BA),2T-\BA)2i] - 0(1). 
Since N-^+^^^^"^^ A21 = 0(1), it is sufficient to prove 

(S^)ii + iBA)i2T-\BA)2i = oil). 

By the LHS is B12A21+ BiiAi2T-'^B22A2i and so we see that it is sufficient 

to show that 

S12 + B11A12T-IS22 = 0(7V-l+l/(2'")). 

Using (j4.32(l this reduces to showing that 

B12 + BiiAi2RC^^RB22 = o(iV-i+i/(2™)) 

or 

RB12R + iRBiiR)iRAi2R)C^iiRB22R) = o(7V-i+i/(2™)). 

Using 7V-i+i/(2™)i?Bi2i? = -iV-i+i/(2m)52i +0(1) which follows as in (|4.36|l . we 
are reduced to proving finally 

(4.37) -B21 + S22v42iCfiiSii = o(A^-i+i/(2™)). 

But 

— -S21 + B22A2iC^i Bii = 

by (taking the transposes of) |DGI (5.12)]. The proof of Proposition H^^ is complete. 

□ 
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Remark 4.1. The second relation in H4.31|l was sharpened recently by Kriecherbauer 
and Vanlessen |KV| who showed that the o(l) terms are in fact identically zero. 
One might hope that this improved result could be used to strengthen the estimates 
in (|4.29() . (|4.30|) . This is indeed the case for (|4.30(l : one can show that 6*12 = £'12 
identically. However we have not been able to use |KV | to improve the estimate in 

4.3.1. The contribution of the correction term to the 12 entry of Kj^^i. In view of 
(|4.27|l . since {SD){x,y) = —dyS{x,y), the correction term has the form 

(4.38) ^lixf ■ Gil • $i(y) + ^lixf ■ G12 • $2(y). 

Again set x = y — r/^-*. Using Proposition 14.21 and proceeding in the same 

way as in Subsection 14 . 2 . 1 1 we find that as iV ^ 00, the term ()4.38|l . multiplied as 
before by (^^^7^373 )^ becomes 

(4.39) _lAi(e)Ai(,7)+o(l)e-^«e-"' 

uniformly for ^,77 > Lq. Note that the sum of all elements of (the binomial matrix 
in the limiting form of) 1^12 is, up to a sign, the same as for 1)21 • 
Remark: Note also that the only new element in the above proof as compared with 
the case /3 = 4 in Subsection 14.2.11 is that the matrix Gn is only asymptotically 
(and not identically) skew symmetric. This leads to the estimate o(l)e~°^e~'^'' in 
place of H4.16|l . 

4.3.2. The contribution of the correction term to the 11 and 22 entries of Kn^i. We 
consider the 11 entry of K^^i (the 22 entry is considered in the same way). Using 
(|1.47|l , H4.27|) we rewrite the correction term as 



(4.40) 



-^i{xf -Gn ■ <Pi{t)dt^ -^,{xf -012- '^2{t)dt 

- <^i{xf ■ Gil ■ e$i(+oo) - <i>i{xf ■ G12 • e$2(+oo). 



Again set x — y — ri^-^\ The first two terms can be treated in the same way 
as in Subsections 14.2.21 and 14.3.11 More precisely we find that the first two terms 
in H4.4(J|) . multiplied by ^^^^2/3 ; become, as iV ^ cx) 



(4.41) -^Ai{0j^Ai{t)dt + o{l)e 

uniformly for ^,77 > Lq. Now consider the (scaled) sum of the last two terms in 



(4.42) 2^i^^^^^iN))T . . gci>i(+oo) + <i>i(^W)^ • G12 • e<i>2(+oo)) . 



By (|4.3|l . (|4.1U|) . Proposition l4.2l this becomes as A^ ^ c» 
(4.43) 



aArAr2/3 ^1/2 (2m)i/2Ari/2 

X {e^ • Gn • (b + 0(1)) + • G12 • (a + o(l))} Ai(0 + o(l) 
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Setting a = e — b, we find that (|4.4,S|I reduces to 
(4.44) 

i Ai(e) + o(l)e-=« + (e^ . Gn • b - • G12 • b) • Ai{0 ■ OliV-i+i/^^")). 
So if we could prove 

(4.45) -Gii-h-e^ ■Gi2-h^o{N^-^^^^"'^), N ^ 00, 
then we would find that H4.44|l equals 

(4.46) iAi(0+o(l)e-^f 

uniformly for ^ > Lo- 

We prove (|4.45|l . We will, perhaps surprisingly, use a property of the (3 = 4 
correlation kernel Sn/2A- is not clear how to prove 14.45(1 directly using the 
asymptotic properties of {D(t)M+j,(t>N+k) and {e(f)N+j,<l>N+k) given in | DG| . More 
precisely, (|4.45|l follows from (|4.29|) , H4.30|l and the relation 

(4.47) b + GiYSiiA2-a = o(l), N ^ 00, 

which is proved by dividing (|1.46(l by (^f^v)^^^ ^^'^ using (|4.:^|l as cx3. Mul- 

tiplying H4.47|l from the left by e-^ • RD21 and noting b = i? • b, a = i? • a, we 
find 

e"^ • RD21R • b + • RD2iC^^BiiDi2R • (e - b) = o(iVi-i/(2m))_ 

But the second matrix is skew symmetric (see (|4.28|l et seq.). By (14.29(1 the above 
relation becomes 

(4.48) e"^ • RD21R • b + • Gil • b = oiN^-^/^"^"^^). 
But from 

i?i?2li? = -i?12+o(A^l-l/(2r„)) 

and hence 1(1301) imply ((05|) . 

Collecting the estimates (14.41(1 . ((4.4611 we see that since J^^Ai{t)dt ~ 1, the 
correction term in the 11 entry has the form 



^-Arnl^l- Aiit)d?j +o(l)e-^« 



(4-49) ' ^ , 

= -Ai(e)/ M{t)dt + o{l)e-^^ 

uniformly for S,,r]> Lq. The correction term in the 22 entry has the same asymp- 
totic form with ^ and rj interchanged. 

4.3.3. The contribution of the correction term to the 21 entry of Kjsi^i. By ((1.49(1 . 
(14.27(1 the correction term in this case has the form 

(^j\i(tfdt^ ■Gii-e<^>i{y)+(^J\iitfdt^ ■ G12 ■ e<^>2{y) 
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which equals 

J\i{tfd?j - Gil •e$i(+oo)+ (^J\i{tfd?j -Gia •e$2(+oo) 



by H1.47I) . Again set x — C*-^-*, y — t?^^-*. A calculation very similar to the one in 
Subsection l4.3.2l using the last estimate in H4.1()|l in place of the first, leads to the 
following asymptotic form for the 21 correction as N —^ oo 

(4.50) J\i{s)ds + ^(^J\i{s)ds^(^J^^ Ai(t)dt^ +o(l)e-" 

uniformly for ^,77 > Lq. (Recall that there is no overall scaling factor for the 21 
entry.) This completes the analysis of the contribution of the correction term to 
the 21 entry of K^^i. 
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